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The breaking of an approximate discrete symmetry, the 
final stages of a first order phase transition, or a post- 
inflationary biased probability distribution for scalar fields are 
possible cosmological scenarios characterized by the presence 
of unstable domain wall networks. Combining analytical and 
numerical techniques, we show that the non-spherical collapse 
of these domains can be a powerful source of gravitational 
waves. We compute their contribution to the stochastic back- 
ground of gravitational radiation and explore their observabil- 
ity by present and future gravitational wave detectors. 
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In order to further our understanding of the physi- 
cal processes that took place early in the history of the 
Universe we should explore potential events which may 
have left an imprint detectable by current or future ex- 
periments. One exciting possibility is that certain pri- 
mordial processes generated a stochastic background of 
gravitational waves. During the past decade or so, several 
cosmological sources (as distinguished from astrophysical 
sources such as coalescing binary systems [0 ) of stochas- 
tic gravitational waves have been proposed H. These 
include inflationary models [p[, cosmic strings |4j , strong 
first order phase transitions fcUfJ , and nontopological soli- 
tons 0. 

All ground-based detectors (interferometric or reso- 
nant bar), such as LIGO or VIRGO, probe the frequency 
interval 10Hz < / < 10 4 Hz. Lower frequencies are to 
be probed in space, as with the planned LISA mission, 
which has a projected sensitivity of 10 _4 Hz < / < 1Hz. 

Here we propose another powerful source of primor- 
dial gravitational waves, the non-spherical collapse of 
bounded domain walls, or bags, which separate differ- 
ent vacuum regions. Zel'dovich et ai, have shown that 
the appearance of domain walls is a direct consequence 
of the breaking of a discrete symmetry H . In the same 
work, domain walls were shown to be incompatible with 
big-bang cosmology, as they would create a power law 
expansion ruled-out by observations. Thus, if domain 
walls were created they had to disappear. One of the 
mechanisms proposed by Zel'dovich et al. was to con- 
sider an approximate discrete symmetry as opposed to 
an exact one; the difference in energy density between 
the two vacua generates a pressure force that eliminates 
the walls. Other scenarios for the disappearance of the 
domain walls were investigated in Ref. [||. 

The existence of walls and other topologically stable 



defects were explored by Kibble, who also suggested the 
possibility of an approximate symmetry as a mechanism 
to eliminate the walls JhJ. Several authors explored 
this idea further [pd|Js[|. Recently, it was shown that 
unstable domain walls could also be created during a 
post-inflationary non-thermal phase transition p"2fl . The 
evolution of the unstable domain-wall network has been 
studied numerically jll| and analytically E3]. 

We can envisage at least three scenarios where unstable 
domain wall networks can be generated: i) the breaking 
of an approximate discrete symmetry in a thermal phase 
transition; ii) biased fluctuations in a non-thermal post- 
inflationary scenario; iii) during the final stages of a first 
order phase transition, as bubbles of the true vacuum 
percolate, leaving a disconnected network of shrinking 
false vacuum domains. 

We can model all these scenarios with a biased double- 
well potential for a real scalar field, 

V(cl ) ) = ^(cf > 2 -cf > l) 2 + hcbo4>( K l^-4) , (1) 

where h > is a dimensionless constant which biases 
the potential towards the +0 O minimum. The energy 
density difference between the two minima is A = ^h<f>Q, 
while the energy density barrier between the maximum 
at 0o = —(h/X)<po and the global minimum (+</>o) is 
B = j<J)q[1 — 0{h/\)]. Coordinates scale as x = ^/Xfo x, 
while energies scale as E = (-s/A/^o) E. 

The formation of a domain wall network during a phase 
transition has been discussed in several works Jl^^,^,|l2| . 
There are two possibilities, determined by the probability 
the field has of landing on either vacuum, p± . There is a 
critical "percolation probability" for a given vacuum, p c . 
For cubic lattices, p c = 0.31, the value we will adopt here 
Jl5t - If both p + and p_ are larger than p c both vacua 
percolate, being separated by a domain wall stretching 
across the lattice. At formation, the wall will be rather 
convoluted, with local average curvature given by the 
typical fluctuation scale, the correlation length £. If only 
P+ > Pc, most of the volume will be in the +- vacuum, 
with isolated clusters (bags) of the negative vacuum with 
a distribution function, f(r) ~ r -1,5 exp[—r], where r is 
the number of cells of unit volume ~ £ 3 in a given 
cluster We will call this case the "no-percolation" 

case. 

The energy density of an isolated domain is given by 

pIM = U^) 2 + lvfa-v<p b + v(0 b ) , (2) 
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where </>&(x, t) is the field configuration describing the bag 
at time t. We define the average radius of an isolated 
domain by P av = (f dV | x | p[fa]) / J dV p[0 b ]. 

Two forces will act on the walls; the tension, Pr, 
which will act to straighten the walls, and the vacuum 
pressure, Pa, which will cause the domains to shrink. 
Here, we will consider two possibilities. First, the "no- 
percolation case" , where the false vacuum bags disap- 
pear very quickly after they form, that is, with Pa > Pr 
right at formation. This case is characterized by having 
several small bags within the horizon. Second, we will 
consider the "percolation case" , where the tension force 
will straighten the walls for a while, up to the horizon 
scale, before the vacuum pressure acts to accelerate the 
walls against each other, causing the formation of large 
unstable bags |J. 

We compute the output in gravitational radiation from 
the collapsing domains within the full, linearized gravity 
approximation More details will be provided in a 

forthcoming publication. In the meantime, the reader 
may consult Ref. , as we follow a similar approach. 

The computation of the spectrum in gravitational 
waves from collapsing 3-dimensional domains in the lin- 
earized gravity approximation is a rather cumbersome 
and computer-intensive task. As the authors of Ref. 
we have limited ourselves here to studying the spectrum 
for situations with axial symmetry. In our case, as toy 
models to more realistic bag configurations, we consid- 
ered two types of domains, ellipsoids and deformed Gaus- 
sians. Contrasting our results with the full 3-dimensional 
calculations for the simpler quadrupole approximation, 
we conclude that here we are providing a lower bound 
for the total output in gravitational radiation. 

The total energy radiated in the direction k into the 
solid angle ft at frequency uj is given by, 

= Geo 2 [T zz (k, w) sin 2 9 + T xx (k, to) cos 2 9 
- T yy (k, to) ~ 2T xz (k, u) sin 9 cos 9} 2 (3) 

where T n (k, lo) are the Fourier-transformed spatial com- 
ponents of the energy momentum tensor for the field con- 
figuration 0b ||. 

From simple scaling arguments, the integrated energy 
radiated in gravitational waves can be written as, 



i?Gw(-Rav) = e(Pav)G(T 2 P: 



2 U 3 
av ' 
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where a is the surface energy density for the domains, 
and e(Pav) = £o ("f^) i s ^ ne efficiency parameter, with 
6q a number to be determined. £o is the T = correlation 
length. 

For ellipsoids, we wrote 



6 fc (r, z) = 4> + + 



(l - tanh(e/v / 2)) , (5) 



where e is the solution of r 2 /(a + e) 2 + z 2 / (b + e) 2 = 1. 
For deformed Gaussians, we wrote, 



6 (r, 9) =(/)+ + 0- - 0+)exp 



[Ro{l + SR)Y 



, (6) 



where SR = ^ z m Pi jm Yz, m (0, tp) measures the distortion 
from spherical symmetry. For axial symmetry we kept 
to = 0. These configurations are then used as initial 
data for solving numerically the Klein-Gordon equation 
for 0(r, z,t). As the configuration evolves in time, we 
compute the emission of gravitational radiation using Eq. 

i 

As shown in Fig. 1, the energy of the initial configura- 
tions can be well fitted by P (P av ) = P° (-^p) • For 

spherically symmetric domains with degenerate vacua 
7 = 2. Writing also E (R, 



= crP 2 v , we obtain that 



the surface tension scales as a 



7-2 



From 



left to right, the sample points shown are: for ellipsoids 
(a b) = [(3 1.5), (2 4), (4 2), (4 6), (6 4), (6 10), (10 6)]; 
for deformed Gaussians (I Rq SR) = 
[(4 3 .5), (3 3 .5), (4 4 .8), (3 5 .3), (4 5 .5), (4 6 .8), (4 8 .5)]. 
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FIG. 1. The initial energy of the nonspherical domains vs. 
their average radius for degenerate (h — 0) and nondegenerate 
(h = 0.15) potentials. The power law fits are shown explicitly. 

In Fig. 2, we show the integrated energy in gravita- 
tional waves, which can be fitted by a simple power law 
relation, [Mpi — 1.2 x 10 19 GeV is the Planck mass] 



^Gw(Pav) — Eq 



gw 



pi 



Pav 
lb" 



(7) 



Comparing Eqs. [| and|^, we obtain that (3 = a+2'y— 1, 
and eo = pv A] . The quantities with a tilde are 

v2(a ) J 

dimensionless quantities and the quantity in the squared 
brackets comes from including thermal corrections to the 
correlation length, £(Tg) = £o/(2\/A) = ^ XTa i evalu- 
ated at the Ginzburg temperature Tq ~ 4</>q/(1 + 4A), 
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the temperature where the wall network forms in a ther- 
mal phase transition UUM ■ 
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FIG. 2. The integrated output in gravitational radiation 
generated during the collapse of the nonspherical domains vs. 
their initial average radius. The sample points shown are the 
same as in Fig. 1. 

Since the time-scale of collapse is well approximated 
by the average size of the domains, the rate of emission 
is, 



M 2 



(3-1 



(8) 



The energy density in gravitational radiation for clus- 
ters with r and r + dr cells in a time interval between t 
and t + dt is 



dpcw — ^gw n{r)dr dt 



(9) 



where r ~ (-Rav/£) 3 i an£ l the number density of r-clusters 
can be written as n(r) — Cr e~ r /Vjt. Demanding 
that the bags occupy a fraction p_ of the horizon vol- 
ume (p_ < p c ), we can fix the normalization constant as 

The typical bag lifetime r is of order R av <C Xh, where 
X H = (45/47T 3 #* y^Mpi/T 2 is the horizon length and 
is the number of relativistic degrees of freedom at tem- 
perature T. The total output in gravitational radiation 
is thus. 
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r /8/3-3/2 
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(10) 



where fp(r max ) = J x max r^/ 3 - 3 / 2 e" r dr, is 0(1) for all 
interesting values of its parameters. 

Using that the Universe is radiation-dominated with 
energy density p la d — (ir 2 / 30) g*T 4 , the fraction in energy 
density from gravitational waves is, neglecting factors of 
0(1), and using T G ~ 20 o , 



J1 GW ~ 10 2 A[2-/X] 




(11) 



Since the domains will mostly have linear dimensions 
of order i? av ~ £, the dominant frequency of emission is 
lu ~ l/i?av — \/2A[2\/A]</>o- Redshifting these quantities 
we obtain, 



h 2 O 



1.7 x 1Q- 
[2VA]" s 



100 
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(12) 



and, f° ~ 1.3 x 10 10 VA[2\/A] (f^J 3 Hz . 

Using that the detector strain at wave band given by 
observational frequency / for stochastic gravitational ra- 
diation is h c (f) = 1.3 x IO-^^Oqw] 1 / 2 ^//) 01, we 
get 



hc(f) 



-30 



4.1 x 10 

[2VA]-V2 



10 4 



100 V 

5* y 



A/pi 



(13) 



The no-percolation case is characterized by a very high 
frequency spectrum and small amplitudes, beyond the 
presently projected sensitivity of ground-based interfer- 
ometers. The situation is quite different for the percolat- 
ing case. 

When both vacua percolate, they are separated by a 
convoluted domain wall, with initial average curvature 
i? av ~ £. The tension force dominates the wall dynamics 
for a while, until the vacuum pressure starts accelerating 
the walls against each other, eventually causing the for- 
mation of large, unstable domains. In order for this to 
happen, the asymmetry must satisfy h < 3 -j^A[2\/A] 3 Sq, 
while, for the walls to disappear as their average curva- 
ture reaches cosmologically significant scales (R av < Xh), 
h > {\) 1 ' 2+ ^- 2 \^/g:(j) lMpif-<E^. To these bounds 
we add the percolation constraint, p_ > 0.31, which be- 
comes, h < 0.15A §. 

Writing the average radius of the domains as R av = 
a\n, the energy density in gravitational waves from col- 
lapsing domains is, (a is not the same as the ellipsoidal 
axis) pew — -%pr, or with k = 45/327T 3 , 



Pgw 



2g° w 4 V- 3 A^- 1 )/ 2 / 0o 



K (3-/3)/2„(/3-3)/2 \M, 
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(14) 



The typical frequency will be, / av ~ [27raA#] 1 , which 
redshifts to 

f° ~ 6.1 x lO 11 ^ 1 ( ^-Y ( -^-) Hz , (15) 

J VlOO/ \MpiJ ' v ; 

while the fraction of the energy density in gravitational 
waves today is 
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(16) 



The detector strain in the wave band / is given by, 

i ?+± 



hc(f) 



a (l-/3)/2 A (l-/3)/4 K -/3/4 



\M Pl 



(17) 



We have summarized our results in Fig. 3, where the 
detector strain vs. frequency are plotted for several val- 
ues of (3 and compared with the LIGO (initial and ad- 
vanced) and LISA sensitivities. The solid squares locate 
an electroweak phase transition at T = 200 GeV. The 
solid triangles, a hypothetical transition at T = 10 9 GeV. 
Continuous curves are for ellipsoids while dashed curves 
are for Gaussians. The dot-dashed curve is for a Gaussian 
with a = 10~ 2 . LIGO is mostly sensitive to collapsing 
walls at (f>o ~ 10 9 GeV or so, while LISA can probe the 
electroweak phase transition for a comfortable range of 
parameters. The spectrum from collapsing bags, being 
very pronounced at a given frequency, is quite different 
from the mostly flat spectra of cosmic strings and in- 
flationary models A more detailed analysis will be 
presented in a forthcoming publication. 
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FIG. 3. Predicted detector strain vs. dominant frequency 
band from collapsing nonspherical domains. We also display 
the sensitivities of the LIGO and LISA detectors. The solid 
cubes are for an electroweak transition at T = 200 GeV and 
the solid triangle for a hypothetical transition at (f>o = 10 9 
GeV. 
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It is clear that the stochastic background of gravita- 
tional waves left behind from collapsing domain walls can 
be detectable for several situations of interest, illustrat- 
ing the potential impact of gravitational wave astronomy 
on applications of particle physics to the early Universe. 
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